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Abstract 


A numerical model has been developed to simulate the effect of combustion zone geometry on the steady state and transient performance 
of a tubular solid oxide fuel cell (SOFC). The model consists of an electrochemical submodel and a thermal submodel. In the electrochemical 
model, a network circuit of a tubular SOFC was adopted to model the dynamics of Nernst potential, ohmic polarization, activation polarization, 
and concentration polarization. The thermal submodel simulated heat transfers by conduction, convention, and radiation between the cell and the 
air feed tube. The developed model was applied to simulate the performance of a tubular solid oxide fuel cell at various operating parameters, 
including distributions of circuits, temperature, and gas concentrations inside the fuel cell. The simulations predicted that increasing the length of 
the combustion zone would lead to an increase of the overall cell tube temperature and a shorter response time for transient performance. Enlarging 
the combustion zone, however, makes only a negligible contribution to electricity output properties, such as output voltage and power. These 
numerical results show that the developed model can reasonably simulate the performance properties of a tubular SOFC and is applicable to cell 


stack design. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


A fuel cell is a device that converts chemical energy into elec- 
tricity directly through the electrochemical reaction of a fuel (for 
example, hydrogen) with oxygen. Compared with traditional 
thermal-electricity generators, a fuel cell is not only more effi- 
cient but also more environmentally friendly because of low or 
zero pollution produced [1,2]. Planar and tubular geometric con- 
figurations are the two typical designs for SOFCs. There have 
been over 40 years of theoretical, laboratory, and modeling stud- 
ies of tubular SOFCs, a longer time than for planar SOFCs [3]. 
The technique of design and manufacture of tubular SOFCs is 
relatively mature [3]. 


* Corresponding author. Tel.: +86 451 8254 2449; fax: +86 451 8254 2449. 
E-mail address: jiajunxi99 @ sohu.com (J. Jia). 


0378-7753/$ — see front matter © 2007 Elsevier B.V. All rights reserved. 
doi:10.1016/j.jpowsour.2007.06.057 


Most experimental investigations have focused on the over- 
all electrochemical performance [4,5], while computational 
modeling has focused on systematic characteristics at vari- 
ous operation conditions. Because of their validity, calibrated 
numerical models are expected to be able to achieve system 
optimization for both steady state and transient performances. 

Various models have been developed to simulate tubular 
SOFC performance [3,4,6—20]. The features of these models 
are different. Some are steady state models [3,4,6—11]; others 
are models [12—20] addressed to transient operations. 

In general, the electrochemical models of tubular SOFCs 
include three polarizations: ohmic polarization, activation 
polarization and concentration polarization. Hirano et al. [4] 
developed a model to examine the roles of ohmic and con- 
centration polarizations on electrochemical properties with the 
assumption that the effects from activation polarization were 
negligibly small. By contrast, the model Aguiar and Chadwick 
developed [7] did not include concentration polarization because 
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Nomenclature 

A area (m) 

C concentration (mol m~?) 

D diffusion coefficient (m? s7!) 

Fact activation energy (J mol`!) 

F Faraday constant (96485 C mol™!) 
AG change in Gibbs free energy (J mol~!) 
i current density (A cm7?) 

ig exchange current density (A m~?) 

i current (A) 

Ji transport rate of specie i (mol m~? s7!) 


heat exchange coefficient (W m2 KJ) 
length of cell (m) 

partial pressure 

cell radius coordinate (m) 

universal gas constant (8.314 J mol—! K7!) 
terminal voltage (V) 

fuel utilization 

oxidant utilization 

electrical power (W) 

cell axial coordinate (cm) 

molar fraction of specie i 

electrons transferred per reaction 


ATR TU 


35 


N > & 


Greek letters 


a transfer coefficient 

E porosity 

n polarization (V) 

À thermal conductivity (W m`! K7!) 
p specific resistivity (Q cm) 
T tortuosity 

Subscripts 

a anode 

act activation polarization 

air air in the air feed tube 

b air feed tube 

c cathode 

con concentration polarization 
CZ combustion zone 

e electrolyte 

8 gas 

ohm ohm polarization 

S cell solid structure 


the principal electrochemical reactions were considered to be 
kinetically controlled. Therefore, only activation polarizations 
were taken into account. For heat transfer, some models have 
simply considered taking heat conduction and convention into 
account to simulate the thermal properties: for example, the 
models developed by Hirano et al. [4] and Campanari and Iora 
[8]. Radiation, which plays an important role in heat transfer 
when the temperature of the cell tube is high [3], has been rarely 
studied. 


Achenbach [21,22] developed a transient model for SOFCs 
that included ohmic, activation and concentration polarizations 
and radiation heat transfer. M.A. Khaleel et al. [23] numeri- 
cally analyzed the performance of SOFCs at steady state and at 
startup using a model that did not take radiation heat transfer 
into account. Both models were used to simulate the operation 
of planar SOFCs, not tubular SOFCs. The simulation results 
from these models showed the overall performance character- 
istics of the modeled fuel cell to be reasonably estimated. The 
computational time was also acceptable. 

Nonetheless, it is clear that little attention has been paid to the 
effects of fuel cell geometry on performance or model design. It 
is also recognized that the effect of the combustion zone geom- 
etry on fuel cell performance has seldom been examined either 
experimentally or numerically. Theoretically, the heat transfer 
through the combustion zone should play an important role on 
fuel cell performance because of the large temperature differ- 
ence produced. For example, when tubular SOFCs work in a 
cell stack, the unreacted fuel can flow into the combustion zone 
to react with the depleted air. The combustion products form the 
exhaust and the heat from the combustion zone can be used to 
preheat the inlet air. 

In this work, a model was developed to simulate the opera- 
tion of a tubular SOFC with a specified combustion zone. The 
electrochemical model includes all three polarizations: ohmic, 
activation and concentration. The radiation between the cell tube 
and the air feed tube was also included in the thermal model. 
The effect of the combustion zone geometry on the steady state 
and transient operation performance of the cell was analyzed in 
detail. 


2. Model description 
2.1. Cell configuration 


Fig. 1 illustrates the configuration of the SOFC stack. A 
tubular SOFC with combustion zone is shown in Fig. 2. 

Fuel enters near the outside of the closed end of the cell 
and flows along the cell to the open end. Air is provided to the 
inside of the cell via an air feed tube and moves to the closed 
end of the cell. Oxygen in the air fed to the cathode accepts 


Fig. 1. Scheme of the SOFC stack. (1) Air; (2) exhaust; (3) fuel; (4) unreact air; 
(5) unreact fuel; (6) combustion zone. 
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Fig. 2. Configuration of tubular SOFC with combustion zone. 


Table 1 

Geometry parameters of the tubular SOFC 

Component Value 
Thickness of cathode 2200 um 
Thickness of anode 100 um 
Thickness of electrolyte 40 wm 


Thickness of air feed tube 1mm 


Outer diameter of cell 22 mm 
Outer diameter of air feed tube 12 mm 
Length of cell 150 cm 


electrons from external circuit to form oxygen ions. The ions are 
conducted through the solid electrolyte to the anode. At the fuel 
electrode, the ions combine with hydrogen in the fuel to form 
water. Electrons flow from the anode through the external circuit 
back to the cathode. Since the electrochemical is exothermic, the 
cell produces heat as well as electricity. 

After exiting the cell, the fuel and air mix and react in the com- 
bustion zone. The combustion products flow along the outside 
of the air feed tube to preheat the incoming air. 

The dimensions of the cell and each cell component are listed 
in Table 1 [24,25]. 


2.2. Electrochemical model 


To simulate the electrochemical characteristics of the cell, 
the tubular cell is divided into slices by the planes perpendicu- 
lar to x-axis as shown in Fig. 2 and denoted by sliced cell. The 
Nernst potential, electrical current, and polarizations of each 
slice are calculated individually. A circuit composed of electro- 
motive forces and cell impedances models a sliced cell. Fig. 3 
shows the cell equivalent electric circuit. 

The cell terminal voltage is constant because the electrode 
materials are good electric conductors and the resistance of Ni 


Nact,c (1) 
Meon, (1) 

E(1) 
Nonm (I) 
Nact,a (1) 
Nona (1) 


Fig. 3. Cell equivalent electric circuit. 


felt used for the electrical contact between cells is sufficiently 
low. However, the driving Nernst potential varies along the cell 
length. The varying Nernst potential produces varying local cur- 
rent densities. The nonuniform current results in locally varying 
voltage polarizations. Each local Nernst potential is reduced to 
the terminal voltage by the sum of the local voltage polarizations. 
Therefore, the cell terminal voltage is given by 


V= Eo- Nact,a — Nact,c — ohm — Ncon,a — Ncon,c (1) 
where Vis the cell potential and Eo is the Nernst potential, which 


is calculated by 


=AG RT 1/2 
Eo = 0 + In PH (Po) 


2F 2F PHO 


(2) 


2.2.1. Activation polarization 

The development of electrochemical reaction requires over- 
coming an activation energy barrier. This phenomenon defined 
as activation polarization can be described by the Bulter-volmer 
equation [26]. 


Pap {exp (=) — exp =£ = eed } 8) 
RT RT 


where a is the transfer coefficient, z the number of electrons 
participating in the electrode reaction, F the Faraday constant, 
and ip is the exchange current density that can be calculated as 


E 
ioa = Ya (z2 2) (=e zo) exp (- sa) (4) 
P0,a P0,a RT 
0.25 
; PO? Eact,c 
= ere exp | — 5 
10,c ye( 2) p ( RT ) (5) 


Values for Ya, Yc, Eact,a» Eact,o could be found from literature 
[26]. 


2.2.2. Ohmic polarization 

Ohmic losses occur because of resistance resulting from the 
flow of ions in the electrolyte and the flow of electrons through 
the electrode. 

Due to the symmetric current flow through the two halves of 
the cell in Fig. 4, it is only necessary to calculate the resistance 


Interconnection 
Electrolyte 


Anode 
—>, 


Fig. 4. Current path of tubular SOFC. 
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Re Re 


Fig. 5. Equivalent electric circuit of the ohmic resistance. 


of a half slice. The total resistance of a cell slice is the parallel 
combination of the two halves. The equivalent electric circuit of 
the ohmic resistance for calculation is given in Fig. 5. 

The cell equivalent ohmic resistance depends on the anode, 
cathode and electrolyte resistances. 

Ohmic polarization is expressed by Ohm’s law: 


ohm = IJ R, (6) 


where R; = ;6;/A; is the ohmic resistance of anode, cathode, 
and electrolyte in each share of the equivalent circuit of Fig. 5, 
A; the respective area of the section where the current flows, 
6; the corresponding current flow length and p; is the material 
resistivity, which is the strong function of temperature as given 
in Table 2 [1,3]. 


2.2.3. Concentration polarization 

The analysis of concentration polarization should begin with 
the analysis of the transport of gases through porous electrodes. 
The electrode concentration overpotential considers the differ- 
ence in gas concentrations between the electrode-electrolyte 
interface and the bulk. Mass transport models inside the porous 
SOFC electrode must be applied to estimate gas concentrations 
at the electrode—electrolyte interface. In this model, both ordi- 
nary and Knudsen diffusions are considered. The model has 
taken into account the physical properties of SOFC material 
like porosity, tortuosity and pore size of the electrode materials. 
The concentration polarization is given by two terms, related to 
the anode and cathode side: 


0 
RT Xy,Xu,0 RT ao 


n Noon a "con c (i ) 
) F X x9 f F 
Ho H20 


Neon = 


where X and X; is the molar fraction of specie i in the bulk and 
at the electrode—electrolyte interface, respectively, which can be 
estimated by the Fick’s law: 


dc; 
Ji = —D;, 
dr 


(8) 


Table 2 
Properties of SOFC components 


Resistivity (Q cm) 


Cathode 0.008114 exp(600/T) 1 
Electrolyte 0.00294 exp(10350/T) — 
Anode 0.00298 exp(—1392/T) 1 
Air feed tube — — 


Pore radius of electrode (um) 


where D;=[t/e(1/Dim+ Viorel ime Dim the mass diffusion 
coefficient of specie i in the mixture, D;, the Knudsen diffusion 
coefficient of species i, £ the porosity and q is tortuosity. 

Applying the Kirchhoff’s law of current, the currents in each 
share of the equivalent circuit can be obtained [27,28]. 


2.3. Thermal model 


2.3.1. Energy balance equations in tubular SOFC 

The cell temperatures influence the electrochemical model 
and these temperatures affect the local driving voltage, polar- 
izations, and heat generation within the cell. A thermal model 
has been developed to determine these temperatures by a finite- 
volume approach. In the electrochemical model, the cell is 
divided into axial sections or slices. In the thermal model, the 
same sections are used and each section is composed of the 
solid structure, air feed tube and flow passages volumes. Energy 
balances equations for such volumes in each section can be 
described as follows. 

The general form of the energy conservation equation for a 
control volume is, 


dEcy _ 
dt 


ÖOcv — Wev + X nihilin = X nihilo (9) 
i i 


where Ecv is the internal energy (J), Qcy the rate of ther- 
mal energy (W) transferred across the control volume, Wey the 
rate of work transferred across the control volume, }>,n;hilin 
the enthalpy gained due to mass flowing into the element, and 
>> nihilout is the enthalpy loss due to mass flowing out of the 
element. For the steady state, left hand term in Eq. (9) is elimi- 
nated. 
The change rate of internal energy is given by: 


a e Cy - (AV) n (10) 
a or ee dt 


where p is the density (kg m~°), Cy the specific heat at con- 
stant volume (J kg! K7—!), AV the elemental volume (m3), T 
the absolute temperature (K), and f is time (s). 

Tracking of the energy flows in the cell proceeds by making 
separate energy balances for the air in the air feed tube, cathode 
gas, fuel gas and solid structure. 

The energy balances for air preheating volume is expressed 


as: 
i gua bg k+1,p, k+l kp, k 
i i ~ Fai +1, +1, , , 
Paircy A Vair —— At — = Son; Ph; PN n; Phi i 
i i 
LK ACEP — TEP) AD 
air4+a b air 
tle Thermal conduction (W m~! K7!) Emissivity 
3/30% 9.6 0.9 
— 27 - 
3/30% 6.23 = 


- 6.04 0.9 
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where Kair is the convective heat transfer coefficient between the 
air preheating and air feed tube and A4, is the heat transfer area, 
p is the time index and At is the change in time from time step 
p totime p+. 

For the full developed laminar flow, the heat transfer coeffi- 
cient is expressed by [29] 
K= wi 

De 

where Ag is the thermal conductivity of the mixed gas, De the 
equivalent diameter. 

The cathode gas energy balance can be written as 


k,p+1 TEP 
c c — łc 
pecy AVe Ai 
k—-1,pņk—1,p k,p,k,p byk, Pp k, 
i i 
Tk yk 


+K Age"? — Tr?) — Zoho; (12) 


where Kpc is the convective heat transfer coefficient between 
the cathode gas and air feed tube and Kg, is the convective heat 
transfer coefficient between solid cell and bulk cathode air flow. 
The last term on the right hand side is the rate of energy accom- 
panying mass transfer of oxygen out of the bulk to the cathode 
solid. 

The time varying energy conservation equation for fuel gas 
can be written as, 


k,p+1 k,p 
i = 
E At 
_ k-1,p,k-1,p kp, kp sok, p k,p 
=X n; Ph -YOn h” + KATP? — TEP) 
i i 


Ik k, k, 
——— (hy? = ho,r) (13) 


where the first and the second term on the right hand side are the 
rates of energy transfer accompanying the fuel mass flow into 
and out of the element, the third term is the convective heat flux 
from the solid cell to the fuel gas, and the last term is the sum 
of the energy accompanying mass transfer of reactant (H2) and 
reaction product (H20). 

The solid cell unsteady energy balance is written as, 


k 1 k, 
T. Pr — T; P 


pscy AV, 


Às cA k Às k k+l 
zr bass P T; P) = z bass i T; P) 


PRIA? = TOP KAT? = TP?) 


Apt Te?y — EY] oh 
Tey + (A/A) 2F 


A Ay, 0,7; — whP 
(14) 


where the former five terms on the right of Eq. (14) are the net 
energy transferred to the cell by heat transfer via conduction, 


convection, and radiation. The sixth term is the net energy addi- 
tion by the convective flux of reacting species to the solid cell 
and subsequent release of their respective reaction enthalpies. 
The last term is the electrical energy generated in the solid cell 
due to electrochemical oxidation of hydrogen. 

In the cell region, the equation for the air feed tube is 
expressed as: 


kK pth ykp 
poch, AV, —— 
At 

= App (TÉTO? — pry ŽP 4 phe _ pete 

= 5, Aap — Fy") — 5 a — Ty ) 


k, k, k, 
+ KairAa(T. PH T, P) + Koc AQT? — T, P) 


Arait = (FY) 
Tey + (Ab/As)(1/e5 — 1) 


(15) 


In the combustion zone, the equations is 


k,ptl „kp 
Poch, A Vp -h 
x At 
= Ad 4 pele phe Ab 4 TEP petle 
= 5, Ane ~ Fy") — 5 Al, — fy ) 


+KaicAd(Tai? — TP) + Koehoal -TP (16) 
where the last term on the right of Eq. (16) is the net energy 
transferred to the air feed tube by heat transfer via con- 
vection between the combustion products and the air feed 
tube. 

In the combustion zone, fuel and air that exit the cell mix, 
combust and are converted to the combustion products flow. 
The mixture temperature is determined from the energy balance 
of the exiting gas. The combustion react occurs instantaneously 
and the combustion heat is generated in a very thin region. There- 
fore, the heat source is added to the first combustion products 
node. 

The time varying energy conservation equation of the first 
node for combustion products can be written as: 


k,p+1 k 
CZ Tez pa = Tez g 
PczCy A Vez At 


= mix, p ; mix, p k,p;k,p k,p k,p 
= Soni hee? — Son; hy’? + KezAc(T,’? — TSP) 
i i 


-nft A HG (17) 


T,mix 


where the first and the second term on the right hand side of 
Eq. (17) is the rats of energy transfer accompanying the exit fuel 
and air mass flow into and out of the element, the third term is 
the convective heat flux from the air feed tube to the products 
gas, and the last term is the combustion heat. 
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The unsteady energy balance of the rest node for combustion 
product is written as, 
k,p+1 k, 
T. P = rap 


CZ 


Pez cy AVez 


At 
k-1,p,k—1,p k,p,k,p 
i i 
k, 
+KerAc(h? — TKP). (18) 


2.4. Boundary conditions 


The boundary conditions for the energy conservation equa- 
tions are as follows. 

As the air feed tube is divided into n sections and the cell tube 
is divided into m sections, accordingly the combustion zone is 
divided into n — m sections. 

At the entrance of the air and fuel, 


TrA = Tair linlet (19) 
TP = Tylintet (20) 


The air preheating becomes the cathode gas at the exit of the air 
feed tube, therefore the continuity between the feed air and the 
cathode gas can be expressed as, 

TR, = TP (21) 


alr 


The outlet end of cell and the inlet of the air feed tube are 
considered adiabatic, that is, 


T” = 7"! (22) 


TË = TPY. (23) 
2.5. Numerical solution algorithm 


In a SOFC, current and temperature distributions are strongly 
coupled. The electrochemical model is solved with a tentative 
temperature profile. The electrochemical model determines the 
Nernst potential, current and electric power. The thermal model 
accepts these results from electrochemical model and calculates 
the temperature of the gasses and solids. There temperatures are 
applied to the electrochemical model for the next calculation of 
cell Nernst potential, current and power. As the simulation pro- 
gressed the model steps back and forth between electrochemical 
and thermal calculations until the convergence is obtained. The 
equations are solved numerically using the Gauss Seidel method. 


3. Results and discussion 
3.1. Steady state characteristics 
In order to investigate the accuracy of the model, the calculate 


V-i curve is compared with the experimental data in reference 
[5] in Fig. 6. 


1 
-—&- Simulation 
0.8 -&- Reference[5] J 
— 0.6 1 
= 
> Ts=1173K 
0.4 Fuel:89%H2+11%H20 J 
i Oxidant : Air 
Uf =0.85 
0.2 Uo=0.25 J 
0 L 1 L L 
0 0.1 0.2 0.3 0.4 0.5 
i(A/cm?) 
Fig. 6. Cell voltage vs. current density. 
Table 3 


Condition and parameters for simulation 


89% H2 + 11% H20 
21% Or +79% N3 


Fuel inlet composition (molar fraction) 
Oxidant inlet composition (molar fraction) 


Fuel inlet temperature (K) 1173 
Oxidant inlet temperature (K) 1173 

Fuel inlet flow rate (mol s~!) 1.96 x 1073 
Oxidant inlet flow rate (mol s~!) 3.73 x 107? 
Operation pressure (atm) 1 


The relative deviation between the simulated voltage and 
experimental voltage in reference [5] is no larger than 5%. Such 
a good agreement for the terminal voltages between the model- 
prediction and the experiment shows that the present model is 
reliable. 

The operation conditions and parameters for simulation in 
the work are listed in Table 3. 

Fig. 7 shows the temperature profile for the fuel flow (Tp), the 
air in the air feed tube (Tair), the cathode gas flow (Te), and the 
combustion product (Tez). The temperatures of air in the feed 
tube and annulus channel increase gradually in the direction of 
flow due to the heat transferred by convection. The temperature 
distribution of fuel is dependent on the temperature of cell tube. 


1400 r r T 


i=0.4A/cm? 


1350F cz=15cm 


1300 
< 
e 
1250 
1200 
11305 25 50 75 100 125 150 165 
L(cm) 


Fig. 7. Temperature profiles of gas in the tubular SOFC. 
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T(K) 


i=0.4A/cm? 
cz=15cm 


0 25 50 75 100 125 150 165 


L(cm) 


Fig. 8. Temperature profiles of solid in cell. 


Due to the exothermic combustion reaction, the temperature 
of the combustion product has a evident increase from the tem- 
perature of mixed gases, then the temperature of it decreases 
gradually because of the convection on the outside of air feed 
tube between the hot products and the feed tube. 

It is important for the SOFC to work under the maximum 
temperature limit, however, it is difficult to measure the temper- 
ature. In the work, the simulated results for the solid temperature 
profiles along the x-axis are shown in Fig. 8. Both of the tem- 
perature of air feed tube and cell solid structure increase and 
then decreases along the x-axis, the two ends of the cell tube 
have a lower temperature than the middle of the cell tube. The 
reduced temperatures at the ends are due to heat transfer to the 
incoming fuel and air. In the combustion zone, the temperature 
distribution of the air feed tube is affected by feed air cooling 
and combustion product heating. 

Fig. 9 shows the local Nernst voltage and current density. 
Both of them decrease along the x-axis. It is because that along 
the stream, the depletion of fuel and oxide has a significant effect 
on the local Nernst voltage. 

The three kinds of polarization losses, which decrease the 
cell potential from the ideal Nernst value to the real value, are 
shown in Figs. 10 and 11. Among the three polarizations, the 
ohmic polarization plays the most significant role in the cathode- 
supported SOFC. The concentration polarization at the anode 


1 

S 0.95 “Ee 
= (5) 
E 3 
5 09 2 
y a 
E z 
g g 
= 0.85f  i=0.4A/cm? 5 

cz=15cm (S) 

0-8 25 50 75 100 125. 180 


L(cm) 


Fig. 9. Nernst voltage and current density profiles. 


—-& ohm polarizaion 


n(V) 


i=0.4A/cm? 
cz=15cm 


0 25 50 75 100 125 150 
L(cm) 


Fig. 10. Ohmic polarization profiles. 


i=0.4A/cm? 
cz=15cm 


L(cm) 


Fig. 11. Activation and concentration polarization profiles with combustion 
zone. 


side is smaller than that at the cathode because the cathode is 
thicker than anode and the reactant diffuse through the relatively 
thick porous cathode to the reaction site at the cathode electrolyte 
interface. The concentration polarization at the anode side is 
negligible with respect to the other losses. 

Fig. 12 shows the molar fractions variation of gases in the 
cell along the fuel path. The consumption of the fuel gas and 


i=0.4A/cm? 
cz=15cm 


molar fraction 


L(cm) 


Fig. 12. Flow chemical composition profile. 


J. Jia et al. / Journal of Power Sources 171 (2007) 696-705 703 


Table 4 

Effect of the length of combustion on the performance of tubular SOFC 
Lez (cm) 0 15 
i(Acm~*) 0.3 0.3 
Ts, max (K) 1260 1272.4 

Ts, min (K) 1246.4 1258.2 
AT=T;, max-T;, min (K) 13.6 14.24 
T;, aver (K) 1256 1268.3 

V (V) 0.7571 0.7602 
W (W) 130.3 130.93 


0 15 0 15 

0.4 0.4 0.5 0.5 
1295.7 1307.4 1341.1 1352 
1277.7 1288.6 1313.7 1323.5 

18.0 18.75 27.4 28.53 
1289.5 1300.9 1331.2 1341.6 

0.6965 0.7000 0.6404 0.6435 

159.9 160.68 183.3 184.62 


oxidant by electrochemical reaction determines the decrease of 
hydrogen and oxygen along the x-axis. Hydrogen and oxygen 
consumption lead to the formation of water vapor. Therefore, 
the molar fraction of water shows an increase along the x-axis. 
The excess air is needed to provide air-cooling in SOFC as the 
allowable temperature rise of the solid cell is limited by the ther- 
mal stress induced in the ceramic cell components. Therefore, 
the change of O2 molar fraction is not so large as that of H2 and 
H20. 

The effect of the length of combustion on the performance of 
the tubular SOFC in given in Table 4. For the fixed average cur- 
rent density, all the other input data are assumed as in Table 3. It 
can be concluded that increasing the length of combustion zone 
leads to a increasing of cell tube temperature. Although the ter- 
minal voltage and cell power also increase due to the increasing 
of the length of combustion zone, the degree of variation can be 
ignored. 

The effects of the length of the combustion zone on the 
temperature of the cell tube and inlet air inside the air 
feed tube at the open end of the cell region are shown in 
Figs. 13 and 14. 

Fig. 14 shows that the longer the combustion zone, the more 
heat is used to heat the inlet air, which lead to the higher the 
temperature of the cell tube in Fig. 13. For the same inlet air 
at the 1173 K, its temperature increase by 0, 18 and 30 K to the 
1173, 1191 and 1203 K inside the air feed tube at the open end 
of the cell on the condition that the length of combustion is 0, 
15, 30 cm, respectively. 
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Fig. 13. Effect of the length of combustion zone on the temperature profile of 


cell tube. 
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Fig. 14. Effect of the length of combustion zone on the temperature of inlet air 
inside the air feed tube at open end of the cell. 


3.2. Transient simulation 


In order to investigate the transient behavior of the tubular 
SOFC with combustion zone, the transient electrical response 
and the temperature profile of the solid structure have been 
simulated as the average current density increase from 0.4 to 
0.5 A cm7? at t=0s as shown in Fig. 15. During the simulation, 
the inlet temperature and flow rates of fuel and oxidant are kept 
constant, all the other parameters are assumed as in Table 3. 

Figs. 16 and 17 show the transient terminal voltage and power 
profile for the cell. The maximum, minimum and average tem- 
peratures of the solid structure are shown in Fig. 18 during 


0.6 T T T T T 


i(A/cmĉô) 


0.4 


0.3 


0 200 400 600 800 900 


t(s) 


Fig. 15. Step increase in cell current density. 
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Fig. 16. Terminal voltage response due to current step increase. 
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Fig. 17. Output power response due to current step increase. 


the transient for i=0.4 A cm7? up to i=0.5 Acm™?. The step 
increase of the average current results in the increase of the 
power output of cell. The temperatures of cell solid structure 
also increase due to the load change. These temperatures show 
a decaying type of growth toward their finial values. The time 
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Fig. 18. Max, min and average temperature response of the solid structure due 
to current step increase. 
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Fig. 19. Terminal voltage response due to current step increase. 
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Fig. 20. Output power response due to current step increase. 


from the old steady state to the new one is about 800s. These 
curves show the same trends as those presented in [21,22] for 
the planar stack. 

For the same step increase in current density from 0.4 to 
0.5Acm~?, the effect of the combustion zone length on the 
transient electrical response is shown in Figs. 19 and 20. For the 
combustion zone which length is 0, 15 and 30cm, the time to 
reach the new steady state is, respectively 1017, 800 and 679 s. 


4. Conclusions 


A model has been developed to simulate the steady state 
and transient characteristics of a tubular SOFC with specified 
combustion zones. The electrochemical model was designed to 
evaluate the ohmic, activation and concentration polarizations. 
The thermal model includes heat transfer by conduction, conven- 
tion, and radiation. The current distributions, gas concentration 
distribution and temperature distribution were also simulated 
by the model. The influence of the combustion zone length 
was analyzed in terms of both the steady state and transient 
characteristics. 

The results showed that ohmic polarization represents the 
highest loss among the three polarizations in a cathode- 
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supported tubular SOFC, at least one order of magnitude greater 
than the other losses. The middle part of the SOFC has a rel- 
atively higher temperature than that at either end of the fuel 
cell. 

It is found that by increasing the combustion zone length, the 
overall temperature of the cell tube increased, while the terminal 
voltage and the output power hardly changed. In terms of the 
transient response, increasing the length of combustion zone 
led to a quicker response of the cell from one steady state to 
another. None of the previously published tubular SOFC models 
have addressed the effects of combustion zone geometry on the 
steady state and transient performance of the cell. 
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